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13 Abstract 

14 There are significant differences in the biology of males and females, ranging from biochemical pathways 

15 to behavioural responses, which are relevant to modern medicine. Broad-sense heritability estimates differ 

16 between the sexes for many common medical disorders, indicating that genetic architecture can be sex- 

17 dependent. Recent genome-wide association studies (GWAS) have successfully identified sex-specific 

18 and sex-biased effects, where in addition to sex-specific effects on gene expression, twenty-two medical 

19 traits have sex-specific or sex-biased loci. Sex-specific genetic architecture of complex traits is also 

20 extensively documented in model organisms using genome-wide linkage or association mapping, and in 

21 gene disruption studies. The evolutionary origins of sex-specific genetic architecture and sexual 

22 dimorphism lie in the fact that males and females share most of their genetic variation yet experience 

23 different selection pressures. At the extreme is sexual antagonism, where selection on an allele acts in 

24 opposite directions between the sexes. Sexual antagonism has been repeatedly identified via a number of 

2 5 experimental methods in a range of different taxa. Although the molecular basis remains to be identified, 

26 mathematical models predict the maintenance of deleterious variants that experience selection in a sex- 

27 dependent manner. There are multiple mechanisms by which sexual antagonism and alleles under sex- 

28 differential selection could contribute toward the genetics of common, complex disorders. The evidence 

29 we review clearly indicates that further research into sex-dependent selection and the sex-specific genetic 

30 architecture of diseases would be rewarding. This would be aided by studies of laboratory and wild animal 

31 populations, and by modelling sex-specific effects in genome-wide association data with joint, gene-by- 

32 sex interaction tests. We predict that even sexually monomorphic diseases may harbour cryptic sex- 

33 specific genetic architecture. Furthermore, empirical evidence suggests that investigating sex-dependent 

34 epistasis may be especially rewarding. Finally, the prevalent nature of sex-specific genetic architecture in 

3 5 disease offers scope for the development of more effective, sex-specific therapies. 
36 
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47 1. Introduction 

48 Sex, and the existence of two sexes, has revolutionised life on Earth. The success of sexual reproduction is 

49 attributed to recombination between parental chromosomes, which accelerates the loss of deleterious 

50 alleles and the proliferation of advantageous ones [1,2]. The difference in gamete size between males and 

5 1 females is a fundamental property of almost all sexual species. Yet sexual dimorphism extends far beyond 

52 this, from cellular and anatomical specialisation to secondary sexual traits such as ornamentation and 

53 behaviour. Furthermore, there are differences in gene co-expression and metabolome networks between 

54 the sexes [3—5]. It is therefore not surprising that in the field of medicine, males and females frequently 
5 5 differ in core features of disease [6]. 

56 The genetic basis of disease has been intensely researched, with the aim of providing improved 

5 7 diagnosis and therapy. Heritable diseases can be classified as being rare with monogenic aetiology (caused 

58 by a single mutation), or common (prevalence 0.1-1%), caused by multiple genetic variants, each with 

59 high population frequency but small individual contribution to disease risk [7,8]. For these genetically 

60 complex diseases and traits, genome-wide association studies (GWAS) have been successful at identifying 

6 1 loci, but the heritability accounted for by main effects, and by polygenic risk score, remains conspicuously 

62 low [9,10]. This deficit is stimulating the consideration of other factors such as the environment and 

63 epistasis [1 1]. Sex differences in the genetic architecture of common diseases have been known for some 

64 time [12], and recent analysis of large GWAS datasets has resulted in an unprecedented rise in the number 

65 of known sex-specific loci for human diseases and quantitative traits (see next section). Whilst this fact 

66 alone should encourage further investigation, evolutionary theory also predicts the existence of sex- 

67 specific genetic architecture for complex traits via sex-specific or sexually antagonistic selection. 

68 Males and females share genes and genetic variation (excepting the Y chromosome), yet often 

69 have divergent optimum conditions for survival and reproduction [13]. A recent model predicts that even a 

70 marginal difference in selection pressure rapidly amplifies the contribution deleterious alleles make to trait 

71 architecture [14]. Additionally, opposing selection pressures on a shared trait creates sexual antagonism, in 

72 which the strength of positive selection for an allele in one sex will allow it to be maintained even if it is 

73 deleterious to the other sex (see section 3) [15,16]. In this review we summarise recent evidence for the 

74 sex-specific genetic architecture of common diseases, and the evidence for sexual antagonism in the light 

75 of evolutionary theory. We also propose new mechanisms by which sexual antagonism can contribute 

76 towards the genetic architecture of disease, and guidelines for the identification of sex-specific genetic 



77 
78 
79 



effects. 



80 



Gilks, Abbot& Morrow [2013) Evolution of sex differences in disease genetics. bioRxiv, http://biorxiv.org/ 



2 



Downloaded from http://biorxiv.org/ on September 18, 2014 



81 2. Evidence for sex-specific genetic architecture 
82 

83 Broad-sense heritability is the proportion of phenotypic variance in a population sample that can be 

84 attributed to genetic variation [17]. With identical genetic architecture, and assuming a common 

85 environment, trait heritability should be equal in male and female samples. However, in a study of twenty 

86 quantitative traits in humans, eleven showed significant sex-bias in heritability [18]. Following a PubMed 

87 literature search, we identified eighteen independent studies in humans (representing thirty -one traits) that 

88 provided separate heritability estimates for males and females, and also stated whether the difference was 

89 statistically significant. A summary of these data is presented in Figure 1, showing that while fifteen traits 

90 did not exhibit significant sex-bias in heritability, thirteen had a higher heritability in females, and three a 

91 higher heritability in males. Although there may be some bias in these studies (non-reporting, non- 
92 independence of traits or prior selection of traits with known sexual dimorphism), they illustrate that the 

93 heritability of complex traits is commonly sex-biased across a range of phenotype classes. 

94 It is well known that non-genetic factors influence differences in heritability, the most obvious 

95 being sex hormones (androgens, oestrogens and progesterones, secreted from the gonads). These can 

96 create systemic differences between males and females for trait expression, which in turn affects disease 

97 risk and heritability, for example the protective effect of oestrogen on heart disease [19]. However, 

98 experiments using hormone treatment and gonadectomy show that some gender differences in phenotypes, 

99 such as immune response, behaviour, and toxin resistance, are not determined by sex hormones but by sex 

100 chromosome dosage [20—22]. This implies that heritability differences are not always caused by sex 

101 hormones, and can be caused by sex-specific differences in genetic architecture, whereby a genetic variant 

102 has a different phenotypic outcome depending on whether it is expressed in a male or female environment. 

103 The molecular genetic evidence for sex-specific genetic architecture is strong. For gene 

104 expression in human cell lines, 15% of SNPs that control gene expression (expression quantitative trait 

105 loci or eQTL) do so in a sex-specific manner, even in the absence of sex hormones [23]. For complex 

106 traits, GWAS have identified many robust sex-specific loci across a range of human phenotypes. These 

107 results are summarised Table 1, which shows thirty -two loci with sex-dependent effects in the twenty-two 

108 traits studied. The majority of the effects were sex-specific (twenty-eight loci significant in one sex only) 

109 although five sex-biased effects were also reported (significant in both sexes but different magnitude of 

110 effect). One opposite effect direction locus has also been reported from a GWAS (for recombination rate 

111 [24]). Model organisms such as laboratory mice Mus musculus and the fruit fly Drosophila melanogaster 

112 have been used successfully in many phenotype-mapping projects, facilitated by the controlled 

113 environment and flexibility of experimental design. Indeed, many sex-specific eQTL have been identified 

114 in the fruit fly D. melanogaster [25] and laboratory mice [26]. For genetically complex traits, sex-specific 

115 loci have been identified for sleep patterns and ageing in D. melanogaster, and for fat mass and skeletal 

116 traits in mice [27-30]. 

117 Gene manipulation studies of model organisms have identified hundreds of genes with sexually 

118 dimorphic effects on disease-related phenotypes. Interestingly, not only have sex-specific and sex-biased 

119 effects been observed, but also sex-reversed and sexually pleiotropic effects (when the same locus affects 

120 different traits in males and females; see section 5). For example, murine vitamin D receptor disruption 
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121 causes weight loss in males but decreased bone density in females [31], and p53 over-expression in D. 

122 melanogaster increases male life-span but reduces that of females [32]. Furthermore, a screen of 1,332 D. 

123 melanogaster P-element insertion lines identified forty-one mutations that had sexually dimorphic effects 

124 on life-span including six that were in opposite directions [33]. Although gene disruption studies do not 

125 precisely reproduce the effect of natural genetic variation, they demonstrate that different pathways can 

126 control the same trait. The question then is: how does this sexual dimorphism in genetic effects arise? 

127 Insights from evolutionary biology are of great value here, since theory about the origin and evolution of 

128 sex differences is well-developed, both on the phenotypic and on the genetic level. 
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130 Figure 1: Comparison of male vs female narrow-sense heritability estimates from human studies. Red and blue- 

131 coloured data points indicate significant differences between the sexes. Grey data points indictate no significant 

132 difference between the sexes. Triangles indicate a binary/qualitative phenotype. Circles indicate a 

133 continuous/quantitative phenotype. la: Drive for thinness, lb Body Dissatisfaction [103]. 2a: Waist diameter. 2b: 

134 Waist-height ratio. 2c: Body-mass index. 2d: Peripheral body fat. 2e: Hip diameter. 2f: Body weight. 2g: Body height 

135 [104]. 3a: Triglyceride serum level. 3b: LDL cholesterol serum level [105]. 4a: Lung FEV1 (forced exit volume). 4b: 

136 Lung D LC o (diffusing capacity). 4c: Lung VC (vital capacity) [106]. 5: Geriatric depression [107]. 6a: Smoking 

137 initiation. 6b: Regular tobacco use [108]. 7: Sleep reactivity (insomnia) [109]. 8: Alcohol dependence [110]. 9: 

138 Subjective well-being [111]. 10: Reading disability [112]. 11: Reading difficulties [113]. 12: Self-esteem [114]. 13a: 

139 Respiratory sinus arythmia. 13b: Heart beat entropy [115]. 14: Tension-type headache [116]. 15: Lower back pain 

140 [117]. 16: Seasonal mood change [118]. 17: Protein C sensitivity [119]. 18a: Drug use. 18b: Tobacco use. 18c: 

141 Alcohol use [120]. 
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Table 1 : Loci with sex-dependent effects on human phenotypes, identified through genome-wide SNP or CNV analysis 



Phenotype 



Sex test 



Separate 
Separate 

Sep & p-diff 
Separate 
Separate 
Separate 

Sep & p-diff 
Separate 
Separate 
Separate 
Separate 

Sep & p-diff 

Sep & p-diff 
Separate 
Separate 

p-diff 
Separate 
p-diff 

Sep & p-diff 

Sep & p-diff 
Separate 
Separate 

Sep & p-diff 
Separate 

GWAMA* 
Separate 

Sep & p-diff 
Separate 
Separate 
Separate 
Separate 
Separate 
Separate 
Separate 
Separate 

Sep & p-diff 



Individuals 
tested 



Gene 



Chromosome band 



SNP 



MAF 



Male effect 1 ' 



Female effect 



Reference 



Mitochondrial DNA levels 
Heart beat rate (QT interval) 
Waist-height ratio 
Waist-height ratio 
Visceral adiposity 
Mitochondrial DNA levels 
Waist-height ratio 

Plasma homocysteine 
Glycine levels 
Crohn's Disease 
Waist-height ratio 
Waist-height ratio 
Recombination rate 

Uric acid concentration 

Sex-hormone binding globulin 

Uric acid concentration 

Waist circumference 

Low-density lipoprotein (LDL) 

Thyroid stimulating hormone 

Naso-pharyngeal cancer 

Waist-height ratio 

Thyroid stimulating hormone 

Pro-insulin levels 

Obesity & Osteoparosis (bivariate) 

Triglyceride levels 

Type II Diabetes 

Type I Diabetes* 

Thyroxine levels (FT4) 

Thyroid stimulating hormone 

Recombination rate 

Height 

Thyroxine levels (FT4) 

Type II Diabetes 

High-density lipoprotein (HDL) 



384 
3761 
175585 
190803 
117857 
384 
175585 
190803 
1679 
3343 
8463 
175585 
175585 
2500 
19578 
28141 
21791 
28141 
199499 
20512 
26420 
1437 
175585 
26420 
27079 
4355 
24273 
149000 
27530 
17498 
26420 
2500 
1625 
17146 
149000 
11528 



MRPL37 
NOS1AP 
LYPLAL 1/SLC30A10 

THNSL2 
RNF144 
GRB14/COBLL1 

CPS1 

ATG16L1 
PPARG 
ADAMTS9 
RNF212 

SLC2A9 
UGT2B15 
ABCG2 
MAP3K1 
HMGCR 
PDE8B 
MICA/HCP5 
VEGFA 
PDE10A 
DDX31 
SOX6 
APOA5/BUD13 
CCND2 
CTSH 
LPCAT2/CAPNS2 
MAF 
1 7q2 1.31 region 

NEDD4L 
NET01/FBX015 
GIPR 
PLTP 



lp32.3 
lq23.3 
lq41 

2pll.2 
2p25.1 
2q24.3 

2q34 

2q37.1 
3p25.2 
3pl4.1 
4pl6.3 

4pl6.1 
4ql3.2 
4q22.1 
5qll.2 
5ql3.3 
5ql3.3 

6p21.33 
6p21.1 
6q27 

9q34.13 
llpl5.1 
llq23.3 

12pl3.32 
15q25.1 
16ql2.2 
16q23.2 

17q21.31 

18q21.31 
18q22.3 

19ql3.32 

20ql3.12 



rsl0888838 
rsl0494366 
rs4846567 
rs2820443 
rsl659258 
rs2 140855 
rsl0195252 
rs6717858 
rsl047891 

rs715 
rs3792106 
rs4684854 
rs6795735 
rsl 1939380 
rsl670533 
rs734553 
rs293428 
rs2231142 
rsl 1743303 
rsl 2654264 
rs6885099 
na (CNV) 
rs6905288 
rs753760 
rs306549 
rs297325 
rs28927680 
rsl 1063069 
rs3825932 
rs6499766 
rs3813582 
rs2668622 
CNP12587 
rs7240777 
rs8 108269 
rs7679 



0.11 
0.29 
0.29 
0.29 
0.35 
0.39 
0.44 
0.44 
0.30 
0.24 
0.40 
0.42 
0.19 
0.33 
0.22 
0.26 
0.30 
0.12 
0.19 
0.38 
0.29 
0.03 
0.45 
0.50 
0.24 
0.20 
0.07 
0.21 
0.30 
0.48 
0.38 
0.20 
0.02 
0.47 
0.31 
0.18 



0.81 
3.08 

ns 

ns 

ns 



ns 
ns 
ns 



+118cM 
-69.4cM 
-0.22 
-0.03 
0.22 
ns 
-4.03 
-0.17 
OR 3-141 
ns 
0.13 
0.04 
OR 1.75 
0.13 
OR 1.08-1.16 
OR 1.13-1.27 
0.02 
0.12 
ns 
ns 
ns 
ns 
ns 



ns 
2.09 
0.06 
0.05 
Z-score 1.5 
0.32 
0.05 
0.05 
0.04 
0.23 
OR 1.48 
0.04 
0.05 

ns 

+88.2cM 
-0.40 

ns 
0.13 
0.03 

ns 
-0.12 

ns 
0.05 
0.08 

ns 

ns 

ns 

ns 

ns 

ns 
0.06 
+124cM 
-8.1% 
-0.08 
OR 1.06-1.14 
1.68 



[121] 
[122] 
[123] 
[124] 
[125] 
[121] 
[123] 
[124] 
[126] 

[3] 
[127] 
[124] 
[123] 
[128] 

[24] 
[129] 
[130] 
[129] 
[124] 
[131] 
[132] 
[133] 
[123] 
[132] 

[88] 
[134] 
[131] 
[135] 
[136] 
[132] 
[132] 
[128] 
[137] 
[132] 
[135] 
[131] 



MAF Minor allele frequency. Value for similar HapMap population sample stated when study sample MAF not available 
t Effect value is for the correlation coefficient p unless otherwise stated. OR Odds ratio, 95% confidence intervals. 
¥ Result of separate-sex analysis of SNPs previously identified in a standard, main-effects analysis. 
* GWAMA 'Genome-wide analysis, meta-analysis' (Magi et al 2010) 
SNP rsl047891 previously known as rs7422339. 



Gilks, Abbot& Morrow (2013) Evolution of sex differences in disease genetics. bioRxiv, http://biorxiv.org/ 



5 



Downloaded from http://biorxiv.org/ on September 18, 2014 



148 3. Sexual antagonism 
149 

150 Sexual antagonism is a genetic conflict resulting from sex-specific selection acting on a shared genome. A 

151 subcategory of sexual antagonism is intra-locus sexual conflict, where a trait is controlled by the same 

152 genes, and is distinct from inter- locus sexual conflict, which concerns reproductive interactions involving 

153 different loci in each sex. Sexual antagonism has now been demonstrated in a wide variety of taxa, 

154 including plants, birds, mammals, and insects [16,34]. Anisogamy (difference in gamete size) is 

155 considered to be the ultimate source of sex-specific selection [35,36], although ecological factors can also 

156 play a role in shaping patterns of sex-specific selection [37]. The fact that males produce many small 

157 gametes and females few large gametes means that reproductive strategies are fundamentally different 

158 across the sexes, which is thought to result in the evolution of sexual dimorphism [38]. However these 

159 divergent phenotypes must be developed from a common pool of genetic information, making it difficult 

160 to simultaneously achieve optimum trait values in both sexes. Thus, for certain traits a conflict will be 

161 maintained and the sexes will be displaced from their optimum phenotypes. For example, when selection 

162 on females was completely removed, experimentally evolved fruit flies became masculinized in a number 

163 of traits demonstrating that males had previously been displaced from their phenotypic optimum by 

164 counter-selection in females (reviewed in [39]). Pedigree analysis of wild animal populations has also 

165 demonstrated a negative intersexual genetic correlation for fitness i.e. genotypes producing successful 

166 males produce unsuccessful females and vice versa [40,41]. 

167 More formally, sexual antagonism occurs when genetically correlated traits have opposite effects 

168 on male and female fitness. In the simplest case, increasing values of a single trait would increase fitness 

169 in one sex and decrease it symmetrically in the other sex (Figure 2a). In this case the trait is assumed to be 

170 positively correlated between the sexes. However more complicated patterns are also possible, such as 

171 opposite fitness effects of different correlated traits (Figure 2b-c) or asymmetric patterns of selection 

172 (Figure 2d). Consistent with this, a recent study demonstrated that human height was likely to be subject 

173 to sexual antagonism: within sibling pairs, men of average height had higher fitness while shorter women 

174 had higher fitness [42]. This means that the fitness effect of a given height-determining allele will be 

175 context-dependent in terms of sex, and that the population as a whole will be unlikely to evolve towards a 

176 shorter phenotype, despite directional selection in females, because of counter-selection in males. One of 

177 the major evolutionary implications of sexual antagonism is the maintenance of genetic variation that is 

178 deleterious to one sex. Although this has not been fully demonstrated at the molecular level, the 

179 population dynamics of a synthetic sexually antagonistic allele in a laboratory D. melanogaster study 

180 accurately follows predictions of mathematical models [43,44]. 
181 
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Sex-specific fitnesses Intersexual genetic 

correlation 




Trait 1 Female Trait 1 



B) 




Female Trait 1 Male Trait 2 Female Trait 1 



C) 




Female Trait 1 Male Trait 2 Female Trait 1 




-^Q2 Trait 1 Female Trait 1 

183 

184 Figure 2: The various forms sexual antagonism can take. Female fitness functions are shown with red lines, male with 

185 blue lines, and the intersexual genetic correlation with black lines. A. The simplest case (also known as intralocus 

186 sexual conflict) is where the same trait has opposite and approximately symmetric fitness effects on males and 

187 females. The intersexual genetic correlation for the traits is high and positive. B. Sexual antagonism can also occur 

188 when different traits have a high positive intersexual genetic correlation, but are selected in opposite directions in 

189 males relative to females. In the unselected sex (broken lines), selection for the trait in question might be weakly 

190 positive, neutral, or even absent if the trait is sex-limited. C. Although no empirical examples of this type have yet 

191 been demonstrated, it is also possible that traits with a strong negative intersexual genetic correlation could be subject 

192 to sexual antagonism, assuming both traits are selected concordantly across the sexes. A negative intersexual genetic 

193 correlation could occur when the same gene product is incorporated in competing alternative pathways. D. It should 

194 also be pointed out that selection pressures need not be completely symmetric. Non-linear relationships are also 

195 possible. 
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196 

197 4. Sexual dimorphism and resolution of sexual antagonism 
198 

199 A sexually antagonistic trait is expected to go through several evolutionary stages (see Figure 3 for more 

200 detail [16,38]). 

201 1. Initially, the trait is monomorphic and under weak stabilizing selection. 

202 2. A change in the physical or social environment causes the previously concordant trait to become subject 

203 to opposite patterns of sex-specific directional selection. This is the most severe stage of antagonism. 

204 3. Sexual dimorphism then evolves, causing the sexes to come closer to their respective phenotypic 

205 optima, but the antagonism is only partially resolved. 

206 4. Finally, the sexes reach their optimal phenotype values and the sexual antagonism is completely 

207 resolved. 
208 

209 Most research to date has focused on sexually dimorphic traits, under the assumption that this dimorphism 

210 is an indicator of sex-specific phenotypic optima. However the stage of the most severe antagonism and 

211 the largest gender load [45,46] is in fact before the trait in question becomes sexually dimorphic, and gene 

212 expression data from D. melanogaster suggested that most sex-biased genes had already reached their 

213 phenotypic optima [47]. In addition, if sexual antagonism results from correlated expression of different 

214 traits across the sexes, monomorphism in a given trait may not be informative about its likelihood of being 

215 subject to sexual antagonism [48]. This speaks in favour of casting a broad net when searching for 

216 sexually antagonistic loci, and not only investigating traits that are already sexually dimorphic. 

217 Proposed mechanisms for the resolution of sexual antagonism include the evolution of sex-linked 

218 modifiers, alternative splicing, or gene duplication [38,49]. Gene duplication is a popular theory as to how 

219 genes can escape sexual antagonism, by allowing each copy to evolve independently for each sex [50]. 

220 Specifically, this would include genes that are activated by sex hormones or have sex-specific 

221 methylation, and are thus expressed at different levels in each sex. There is debate about the time-scale of 

222 the resolution of sexual antagonism [49,51—55], but regardless of whether the process is fast or slow in 

223 evolutionary time, the outcome is always sex-specific genetic architecture. In this sense, sex-specific 

224 genetic architecture in disease is likely to be an indirect result of past sex-specific or sexually antagonistic 

225 selection. Sexual antagonism or sex-specificity could contribute towards disease risk directly, in which 

226 case the sex-specific or sexually antagonistic alleles are themselves a risk factor [14], or indirectly, 

227 because in its resolved state, sexual dimorphism creates its own intrinsic sex-specific risk factors, such as 

228 behavioural and anatomical specialisation. 
229 
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Trait 1 Trait 1 



230 

231 
232 

233 Figure 3: Predicted stages in the resolution of sexual antagonism. A. Initially, the trait is monomorphic 

234 and under weak stabilizing selection. B. A change in the physical or social environment causes the 

235 previously concordant trait to become subject to opposite patterns of sex-specific directional selection. C. 

236 Sexual dimorphism then evolves, causing the sexes to come closer to their respective phenotypic optima, 

237 but some antagonism remains. D. Finally, the sexes achieve their optima and the antagonism is completely 

238 resolved. Redrawn after information presented in [16]. 
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239 5. How sex-specific selection affects disease architecture 
240 

241 One implication of sexual antagonism is the maintenance of deleterious genetic variation at higher 

242 population frequency than would be expected from mutation-selection balance [43,44]. This leads us to 

243 consider its role in susceptibility to common, genetically complex disorders. Consistent with this 

244 reasoning, mathematical simulation predicts that alleles that are under sex-differential selection (including 

245 sexually antagonistic ones) contribute disproportionately to genetic variation underlying disease 

246 phenotypes [14]. We now discuss in greater depth how sexual antagonism for standing genetic variation 

247 might contribute to the genetic architecture of complex traits. 
248 

249 Unequal endophenotype outcome 

250 For common, complex diseases, risk loci are unlikely to cause disease directly, rather they affect a 

251 quantitative trait that is biologically linked to the disease (an endophenotype) and becomes a risk factor 

252 when the trait value exceeds a certain threshold [56]. This is derived from the endophenotype hypothesis 

253 of psychiatric disorders [57], but can be extended to many disorders, examples of which include the 

254 relationship between adiponectin level and Type 2 Diabetes [58], and triglyceride level with coronary 

255 artery disease [5 9] . 

256 For a correlated trait (with the same genetic architecture between the sexes) an extreme 

257 endophenotype value might be beneficial to one sex but detrimental to the other because of other factors, 

258 such as sexual dimorphism in physiology, behaviour or environmental exposure. In practice, the extreme 

259 trait value might not be beneficial to either sex, but as long as it remains neutral or weakly deleterious, 

260 then the causative alleles persist. One example could be for dyslipidaemia that increases risk of 

261 myocardial infarction in men but not in women, likely due to the anti-oxidant effects of oestrogen [60]. 
262 

263 Equal disease risk but with unequal fitness effects 

264 Fitness comprises both survival and reproductive components. One might implicitly assume that the 

265 reduction in fitness caused by disease is due to disease-related reduction in survival. However, the effects 

266 of disease on the second major component of fitness, reproduction, are not always equal between the 

267 sexes. One example of this is schizophrenia, where males have a consistently greater reduction in 

268 reproductive success than females [61—63]. A second example is for congenital hypothyroidism, 

269 associated with loss in fecundity in women but not in men [64]. These examples illustrate how although 

270 the genetic architecture of disease may be the same, the fitness effect on each sex as a result of the disease 

271 differs. 
272 

273 Sex-specific migration 

274 It has been suggested that the genetic variation for a sexually antagonistic trait may vary between 

275 populations [44], and thus immigration results in the introduction of novel, sexually antagonistic alleles 

276 into the host population. Sex-biased immigration will cause alleles that are beneficial to that sex (and thus 

277 under net positive selection) to be rapidly introduced into a host population, only for the opposite sex to 

278 inherit novel deleterious alleles. 
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279 The same principle could be applied to resolved antagonism. For example, methylation is both 

280 sex-specific [65,66] and population-specific [67]. It is also proposed as a means by which sexual 

281 antagonism can be resolved because it prevents a deleterious allele from being expressed in one sex. The 

282 sex-specific migration results in novel combinations of methylated genes increasing the prevalence of 

283 extreme (deleterious) phenotypes. 

284 Although obtaining empirical evidence for these processes may be challenging, there is good 

285 evidence for large-scale, sex-specific migrations amongst historical human populations from Central Asia 

286 [66,68], the Iberian Penninsula [69], the British Isles [70,71], Central Africa [72], Indonesia [73], and 

287 globally [74,75]. Furthermore, these mechanisms could provide a novel explanation for the outbreeding 

288 depression observed in wild animal populations. 
289 

290 Sexually antagonistic pleiotropy 

291 We define sexually antagonistic pleiotropy as the deleterious effect of an allele on a fitness-related 

292 phenotype in one sex, with a gain in fitness in the other sex through a different phenotype (Figure 2b-c). 

293 One example of this comes from quantitative genetics in which cholesterol levels in men are inversely 

294 correlated with height in women [76]. Body height in humans is sexually antagonistic, with high values 

295 increasing male fitness but reducing that of females [42]. Thus, selection for shorter females causes a 

296 maladaptive response in males by raising their cholesterol levels. Indirect empirical evidence indicates 

297 that pleiotropic genes are indeed less able to escape sexual antagonism [47,77], and thus the involvement 

298 of pleiotropic genes in disease risk seems likely to be amplified by sex-specific selection. 
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299 6. Methods for identifying sex-specific genetic architecture 
300 

301 For SNP -based association testing, the basic approach to identifying sex-specific effects is to analyse each 

302 sex separately, i.e. sex-stratified. In comparison to joint tests, this approach is limited due to the loss in 

303 power caused by partitioning of the sample [78]. A common follow-up to the sex-stratified tests is to 

304 determine whether the association statistics for each sex are significantly different from one another. 

305 Many main-effect studies incorporate sex as a covariate into the analysis, i.e. they are controlling for the 

306 effect. However, whilst this approach acknowledges sex-effects it doesn't allow for their detection. For 

307 binary traits with a prevalence of less than 1% inclusion of known covariates actually reduces power [79]. 

308 A joint analysis that incorporates a genotype-by-sex interaction term tests the difference in 

309 allele frequencies between male and female cases, given their allele frequencies in controls. It is thus more 

310 suited to identifying genetic differences in trait architecture between males and females rather than for 

311 main effects. The regression model with which to test for genotype-by-sex interactions in an unrelated 

312 population sample, is: Y[G,S]=Po+PgG+PsS+Pg x s(G x S)+£, where Y is the phenotype value, G is the 

313 genotype, S is the sex, p is the standardised regression coefficient of each variable, and £ is the error [80]. 

314 Other covariates, such as those used to correct for population stratification, can also be incorporated into 

315 this model. The tests can be performed using open-source software, e.g. PLINK [81] and GenAbel [82], 

316 although to the best of our knowledge only one GWAS [83] and two (candidate gene) studies have done so 

317 [84,85]. For family trio data, a joint, interaction analysis is also possible, exemplified by use of a custom- 

318 designed case/pseudo-control test that detected two loci for autism risk [86]. Meta-analysis of GWAS data 

319 is a powerful and routine approach to increase power in large heterogeneous sample collections, and an 

320 algorithm has been developed in which both sex-specific and main effects can be tested for in a meta- 

321 analysis without loss in power [87,88]. An alternative approach (developed for GxE or GxG interactions 

322 but applicable to GxS) is to use only loci for which the variance differs significantly between genotypes 

323 [89]. This reduces the number of tests whilst enriching for loci most likely to have an interaction 

324 component. 

325 The statistical behaviour of genotype-by-sex tests must be assumed to be similar to genotype- 

326 by-environment tests, in which the interaction term is binary, obligatory and ideally has the same 

327 distribution in case and control populations. Power calculations can potentially be undertaken using such 

328 software as Gene-Environment iNteraction Simulator (GENS) [90] and GxEscan [91]. Analytical hazards 

329 when using an interaction term include artifactual population substructure [92] and incorrect control of 

330 confounders [93] other than sex, such as age, ethnicity, or socio-economic background. Interesting 

331 opposite effect direction effects may arise [24,94] but are hard to interpret without replication or biological 

332 validation. 

333 As more sex-specific analyses of GWAS datasets are performed, it would be informative for 

334 authors to present sex-specific values for i) Trait heritability, ii) The phenotypic variance accounted for by 

335 significant SNPs, iii) Genomic prediction/Risk profile score. Finally, given the extent of sexually 

336 dimorphic interaction networks [4,95,96], pathway enrichment and epistasis testing should be rewarding. 
337 

338 
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339 7. Conclusions 
340 

341 Despite sharing genetic variation, there are profound biological differences between males and females. 

342 This can result in different optima for shared traits, sexual antagonism, and sexual dimorphism. Sex- 

343 specific selection on an allele can have important effects on its maintenance within a population, 

344 potentially allowing deleterious, disease-associated alleles to persist [14,43,44]. This predicts the existence 

345 of sex-specific architecture, and indeed the recent analyses of large GWAS data sets has brought about an 

346 unprecedented rise in the number of robust sex-specific effects in traits of medical relevance (thirty-four 

347 loci for twenty-two traits). In fact, we are aware of only one 'sex-sensitive' GWAS that did not reach 

348 genome-wide significance [83]. Sex differences in disease presentation are often stated as the reason for 

349 investigating sex-specific genetic effects, but given that sexual dimorphism is a resolution of sexual 

350 antagonism, sexually monomorphic traits are more likely to harbour unresolved conflict and thus also 

351 have sex-specific genetic architecture. 

352 Although we have partially excluded sex hormones from sex differences in genetic architecture, 

353 they are a major driver of sex-specific gene expression [4]. Furthermore, as their mode of action is to 

354 activate transcription factors, SNPs that alter binding-sites for sex hormone-induced transcription factors 

355 will have a sex-dependent effect on gene expression. This mechanism is likely to explain some of the sex- 

356 specific GWAS results. Alternatively sex-specific epigenetic modification e.g. methylation [97], will 

357 inhibit gene expression, masking any functional variation in one sex but not the other. One example of this 

358 is known for the ZPBP2 gene and asthma [98]. 

359 Much of what is known about sexual antagonism has been obtained through studies on wild and 

360 laboratory animal populations, as well as mathematical modelling. Identification of the molecular genetic 

361 basis of fitness and of sexual antagonism in model organisms would not only confirm the empirical 

362 observations but also provide a grounding for studies of sex-specific genetic architecture in humans. 

363 Equally so, ecological studies in humans could also provide interesting perspectives, for example how 

364 ecological factors influence selection on specific traits to produce varying degrees of sexually concordant 

365 or sexually antagonistic selection across populations [99]. There is evidence from divergent species of 

366 weak sex-specific trans-eQTLs [26,100], sex-specific residual genetic architecture [101] and sex-specific 

367 epistasis [33]. These studies indicate that many modifier loci for common, complex disorders could be 

368 sex-specific. One example of this is age-at-onset of Parkinson's disease being reduced in males only, by 

369 the catechol-O-methyltranferase gene Vall58 allele [102]. Many monogenic disorders originate from 

370 mutations in sex-linked or mitochondrial genes that, because of their transmission dynamics, are under 

371 sex-dependent selection. However, the role of common genetic variation on these chromosomes in 

372 complex traits is limited by lack of coverage of genotyping chips and suitable analytical methods. 

373 We anticipate that analysis of GWAS data with respect to sex, encouraged by both evolutionary 

374 genetics and recent results presented in this review, will generate many more significant findings and 

375 highlight the potential role of sex-specific and sexually antagonistic selection as a potent force in human 

376 genetic architecture. Finally, we hope that the identification of sex-specific genetic aetiologies in what 

377 otherwise appears to be the same disease will result in the development of more effective, sex-specific 

378 therapies. 
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